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We analyse fall-apart decays of polyquark (tetra, penta and molecule type) hadrons within the 
constituent quark picture. For processes in which a polyquark hadron goes to final states containing 
a light pseudoscalar meson the constraints given by chiral symmetry are implemented. As an 
application of the approach developed, fall-apart decays of a(980) and X(3872) are studied, assuming 
these mesons are polyquark hadrons. Two extreme options - confined diquark-diquark states and 
molecular states - are considered. For a (980), the observed width can be obtained assuming that 

\Q , this meson is a diquark-diquark composite with a relatively large size of around 1-^1.5 fm. The pure 

KK molecular-type state, however, can be excluded. For the X(3872), a sufficiently small width 
can be obtained if it is a dominantly isospin-0 diquark-diquark composite with a very large size of 

fNj , > 2.5 fm. The pure molecular option appears possible if the binding energy is tiny, Eb < 0.2 MeV, 

£jjy corresponding to a huge size. 
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We define polyquark hadrons to be tetra, penta and molecule type hadrons which can be viewed as composites of 
, massive constituent quarks together with antiquarks. If occuring at all, such states can only exist at low excitations 
and with masses close to the sum of the masses of their constituents. The question of the possible existence of such 
states have recently received much attention. In spite of the great activity in this sector y|,|2L|j] many experimental 
and theoretical issues still remain open. Let us recall some of them. The theoretical understanding of the light 
' scalar mesons a (980) and /°(980) is still contradictory: Following earlier suggestions £| the authors of Ref. [{| 
give arguments in favour of the diquark-antidiquark picture of these states. However, the authors of Ref. find 
agreement with the data by using the ordinary qq composition of these scalar mesons if taken together with only a small 
four-quark admixture in form of a loosely bound KK component. The newly discovered heavy meson X(3872) has 
properties which make it unlikely to be an exited charmonium cc state. Instead, it could be a diquark-antidiquark 
I system or may have an important four-quark component in the form of a DD* molecule Q. Also, several newly 
found states with open charm Q may find their explanation as admixtures of usual hadronic states with four-quark 
Q_i' composites @- On the other hand, the existence of a five-quark exotic composite at 1530 MeV, the pentaquark 6, 
appears now less probable according to the negative results of new high statistic experiments . 
' Polyquark hadrons with a composition of four or more constituent quarks are worthwhile to study even if they 
, appear only as components of otherwise conventional hadron resonances: These states or component of states have 
an interesting structure and their hadronic decays proceed by a fall-apart mechanism. A characteristic feature of fall- 
apart processes is that the number of constituent quarks contained in the initial hadron is equal to the total number 
■ of constituents in the final hadrons. The decay proceeds by a rearrangement of the quarks in the initial state. For 
instance, a quark and an antiquark from different clusters composing the initial polyquark state can combine to form 
a meson which then leaves the interaction region. This is quite different from the decays of usual hadrons, in which 
at least one new pair of light quarks is generated. One can expect that the amplitude of fall-apart processes depend 
strongly on structural details of the polyquark hadron, in particular, on the spatial distribution of the constituent 
quarks. 

In this article we will set up the formulas for the fall-apart decay amplitudes. These are approximate equations 
because of the approximate nature of the concept of constituent quarks. The approximation allows a convenient Fock- 
space representation of the hadrons involved where all soft gluon effects and the effects of virtual meson exchanges 
can be viewed as being incorporated into the masses and wave functions of the constituent quarks. We discuss decays 
of polyquarks to light pseudoscalar and vector mesons. Chiral symmetry and the connection of vector currents with 
the vector mesons allow to reduce the decay amplitudes to current matrix elements between the polyquark particle 
and the final hadron which has two quarks less than the polyquark. The dependence of the decay widths on the 
form and size of the quark distributions inside the polyquarks provides useful insights. Following our previous work 
[To| , we first assume that diquarks play an essential role for the structure of polyquark states [Til Il2| . They are 
known to be important for low energy processes 01 . According to the results of |13|. the size of the diquarks made 
of light quarks is taken to be close to the pion size. As in our previous analysis of the hypothetical pentaquark [To| . 
we show that the decay widths of polyquark states are generally quite large. However, the widths can be suppressed 
by assuming a large spatial separation between the diquarks, which then corresponds to an almost molecular type 
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structure. Alternatively, we consider another extreme, namely polyquark hadrons as bound systems of conventional 
hadrons, i.e. truly molecule like states. Also in this case it is found, that a small decay width can only be obtained 
for a relatively large separation of the constituents equivalent to a small binding energy. 



2. THE EMISSION OF PSEUDOSCALAR MESONS AND THE AXIAL CURRENT 



For processes containing a light pseudoscalar meson in the final state one can take advantage of chiral symmetry 
by using the divergence of the almost conserved axial vector current as an interpolating field for this meson. The 
problem of calculating decay amplitudes then reduces to the calculation of matrix elements of this current between 
an initial and a final hadron. 



A. The axial current for constituent quarks 



In QCD hadrons consist of a multitude of almost massless quarks and gluons. However, in the region of low energy, 
hadrons can be described as consisting of only few but massive constituent quarks. Hypothetically, we take this 
simplified picture, with sea quarks and soft gluons considered to be integrated out, also for the polyquark hadrons. 
Because of chiral symmetry, the light mesons - the Goldstone particles of this symmetry - play a dual role: Since they 
have small masses they are tightly bound states of constituent quarks but can, just because of this strong binding, be 
counted as additional degrees of freedom |l4j . The axial vector current obtained from the corresponding Lagrangian 
for consituent quarks and light mesons is well known (see e.g. Eqs. (12), (13) of Refs. [3). To perform actual 
calculations of matrix elements of this current using the constituent quark picture for initial and final states one still 
has to go one step further: For the hadron matrix elements we are considering, the effects of virtual meson exchanges 
should be incorporated into the wave functions of the constituent quarks. The mesonic part of this axial current 
must then be replaced by the constituent quark field operators with the help of the equation of motion of the meson 
fields. In doing this, we will stick to a bilinear form in the constituent quark fields only. The axial vector current in 
terms of constituent fields is easily seen to be a nonlocal one. We write it in form of an operator equation with the 
understanding that it is to be sandwiched between hadron states formed of constituent quarks only. For the example 
of a strangeness-changing transition u — ► s, the axial vector current operator reads 

Sw 5 U - (m s + mu) Q + S 75 U j . (2.1) 

Here S and U are constituent quark fields with constituent masses rag and mu, respectively. ~ 0.7 -r- 1 fl6| is the 
axial coupling relevant for constituent quarks. The divergence equation for the first part of this axial current obeys 

8^ (S^ l5 U) = (m s + mu)Si 75 U. (2.2) 

It holds in the subspace of hadrons formed from constituent quarks, because the effective interaction part of the 
Lagrangian for these quarks should still be chirally symmetric. In the chiral limit (m^ — ► 0) the axial vector current 
(|2.1|) is conserved by virtue of l|2.2|l . 

The divergence of provides an expression for the interpolating field of the pseudoscalar meson: 

*K= 9 Q A mU + mS n ' 2 Sv»U. (2.3) 

The amplitude for a decay with the emission of a if meson is therefore obtained by sandwiching the constituent quark 
field operator 

Q m v + m s = 

t = 9a j SvysU (2.4) 

JK 

between initial and final Fock states formed by constituent quarks. 

In usual form factor calculations (e.g. for quasi elastic transitions) the current operator annihilates a quark and 
creates another. In fall-apart transitions, however, the same current annihilates a quark and an antiquark. For the 
latter process it is therefore convenient to reexpress one quark field by its charge conjugate field. In our example with 
constituent up and strange quarks we obtain 

S^ l5 U = iS T j°j 2 ^-f 5 U, (2.5) 
S l5 U = i5 T 7 7 2 75 C/. 
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Here S denotes the charge conjugate field S = CS T with C — ij 2 j°. 

Constituent quarks do not move fast inside hadrons. Since we will always work in the rest system of the polyquark 
and since the transition operator acts exclusively on the polyquark hadron one can use nonrelativistic expressions to 
rewrite (|2.5[) provided the velocity of the final hadron in the current matrix element is also not large. For simplicity 
we take for the nonrelativistic two component fields the same particle name as in the rclativistic version and denote by 
<j k the two by two Pauli matrices. The current components for conventional and fall-apart processes are now different. 
We denote them by and J^, respectively. For the conventional case one has 

J a = - -^—d k S^a k U : (2.6) 

2vrri(j 2ims 

J A = S ] aU, 

h = T7^—^a k d k U + -^—dkS^akU. 

limu Aims 

For the fall- apart processes the following structures are relevant: 

jo = -iS T a 2 U, (2.7) 

J A = --^—S T a 2 aa k d k U - —— d k S T a 2 a k a U, 
Irajj 2ms 

J5 = iS T a 2 U. 

In the expressions for J° and J5 the products of the lower components of Dirac spinors of order d 2 /tuq are neglected 

in this nonrelativistic approach for the transition operator. To this accuracy we have J\ = — J5. With the help of 
these formulas we can now express all current matrix elements and thus all axial form factors in terms of matrix 
elements calculable in nonrelativistic constituent quark models. We note that the matrix elements of J a, J a an d J5 
do not involve small Dirac components. In the Fock space representation of the hadron states they are obtained from 
the overlap of the wave functions with no derivatives. 

It is seen that the fall-apart transition amplitudes for the emission of pseudoscalar mesons can simply be calculated 
from H2.4fl using J5. Nevertheless, because the calculation involves the constituent quark masses and model hadron 
wave functions, also other form factors should be calculated in a given model. They can provide a consistency check 
for the analysis. As an immediate consequence of the divergence equation (12.21) and Eq. 1)2. 7J1 one finds a constraint 
on the masses of the constituent quarks: The effective quark masses have to match the difference beween the mass of 
the initial polyquark hadron and the mass of the final hadron 

Mi - M f ~ mu + m s . (2.8) 

This follows by observing that in the center of mass system one has for the energy transfer q° ~ M, — Mf, and the 
spatial divergence for slowly moving constituent quarks is small in comparision. 



B. Fall-apart amplitudes for scalar polyquark mesons 

Let us consider the Mi(0 + ) — > Mf(0~) transitions with the emission of a K + meson induced by the strangeness- 
changing axial vector current sj^^^u. A possible application could be the decay of the a + (980) to K + K° and to 
7?7r in case the dominant part of the a + is a (us)(ds) state formed of two scalar diquarks or a molecule or cusp type 
(us)(ds) state formed of two K mesons. 

We start by defining the form factors of the axial vector current 

(M f (p')\sr^u\M. l (p)) - ,9i (q 2 ) (p + p'r+ 9 2(q 2 ) (2.9) 
q = p — p'. Since the K meson pole occurs in the form factor g 2 we define the residuum function r(q 2 ) by setting 

rf) = 2 • (2-10) 

— q z + m K 

The Mi — » MfK + decay amplitude can then be expressed in terms of r(m%) ^3 

A(Mi -> MfK) = - % r -^p±. (2.11) 

JK 
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By rewriting now the strangeness-changing axial vector current sj^^u, by virtue of Eq. I|2.1[l . in terms of the 
constituent quark field operators S and U, we obtain 

(Mfip'^sn^hip)) = g^Mfip'^Sj^UlMiip)) - (m s + ma) f_ 2 ^(M/frOi^s Wfr)). (2.12) 

q + m K 

The matrix elements on the right-hand side can be expressed in terms of invariant functions: 

(M / (j/)|S 7 /i TlsW(p)) = G 1 (q 2 )(p + p'r + G 2 (q 2 )q^ (2.13) 
(M f tf)\S<y 6 U\Mi(p)) = G 5 (q 2 ). 

In the g 2 -region of interest these form factors Gi are now regular functions without poles. The connections between 
g's and G's are 

9l (q 2 ) = g^G^q 2 ), (2.14) 

= 9a [ G m ) — — r G 5W ) 



— q z + m K \ -q~ + m K 

At the pole one finds 

r(m 2 K ) = -g% {m s + m v ) G 5 {m 2 K ). (2.15) 
The divergence equation for the first part of the axial-vector current for constituent quarks (|2.2|l gives 

G 1 (q 2 )(M 2 - Mf) + q 2 G 2 (q 2 ) = -(m S + m^G.iq 2 ). (2.16) 

At q 2 = m 2 K we have 

r(m 2 K ) = g Q A (G 1 (q 2 )(M 2 - M 2 ) + m 2 K G 2 {m 2 K )) . (2.17) 
Gi , G 2 and G5 can be calculated from specific components of the left-hand side of (|2.13() : 



B° = (Mfi-filS^UlMi), (2.18) 
B L = ^(Mfi-q^lS^Um, 
B 5 = (M f (-q)\Sj 5 U\Mi). 



One finds 



For the decay amplitude one gets 



Gi = i( B °-fH' (2 - i9) 

G2 = m { Bo + (2M m7 )bl 

G5 = B 5 . 



9 i { ^P^iG 5 , (2.20) 

JK 



in accordance with 12.11|l . Alternatively, the decay amplitude can also be obtained taking r(m 2 K ) from l|2.17|l and 
calculating G\ and G 2 from l|2.19|l . As long as the divergence equation 12.16fl is respected in our model, the equivalence 
of Eqs. (|2.15|l and l|2.17|l is evident by taking q 2 G 2 from (|2.1til) . The divergence equation itself expressed in terms of 
G±, G 2 and G5 requires for its validity m 2 K = [mu + ms) 2 = (Mi — Mf) 2 together with the nonrelativistic relation 
Bo = —B 5 (in the rest system of the tetraquark). 
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C. Fall-apart decay of the scalar tetraquark a(980) 

As an application of the above formalism, we discuss the decays a(980) — > ttij and a(980) — * KK. We consider only 
two extreme options for the composition of this particle. 

A: a(980) is a confined composite system of two spin-zero diquarks in an iS-state. Then the a + meson has the 
structure 

a+ = (S T ia 2 U)(S T ia 2 D). (2.21) 

B: a(980) is a weakly bound 5-state of two if-mesons: 

a+ = (S T ia 2 U)(S T ia 2 D). (2.22) 

The transition amplitude is obtained by calculating the form factor G5 defined in l|2.13[l using the expression for 
the pseudoscalar current as given in H2.7|l . Changing now to a nonrelativistic normalization of the state vectors, we 
introduce the dimensionless form factor g a +^p(q )'■ 

g a +-,p(q 2 ) = 9 Q a {P{-q\Q T i<J 2 U\a+{p = 0)), (2.23) 

where Q = D, P = rj s for the a+(980) -> tt+tj decay, and Q = S, P = K° for the a+(980) -> K+K a decay. This 
form factor determines the a + decay amplitudes and the corresponding decay rates. The Fock-space representations 
for the tetraquark and the final meson states, as well as the formulas for the corresponding transition amplitudes and 
decay rates, are given in Appendix A. Since we work within a nonrelativistic approach, we set </? = 1. 

For numerical estimates we parametrize the radial wave functions of K, 77, n, and the diquarks by a simple Gaussian 
ip(r) ~ exp(— r 2 /2a 2 ) with the size parameters a„ = ao = 0.9 fm [13, c^k/o-tt = 0.9, a ns /a n = 0.8 ^8|. These 
parameters lead to a good description of the elastic form factors of pseudoscalar mesons at small momentum transfers. 
Because the nonrelativistic approach is used, these wave functions do not provide correct values for the decay constants 
of pseudoscalar mesons. For the decay constants which appear in the interpolating currents (|2.3(l . we therefore use 
their empirical values. We believe this is the right way to proceed since the key quantity calculated in our approach 
is the form factor giq 2 )- 1 

For the case A in which the a(980) meson is a diquark composite we use a Gaussian form also for the motion of 
the two diquarks: 

Vv(r)~exp(-r 2 /2a 2 ). (2.24) 

The corresponding size parameter a a determines the mean distance between the center-of-mass positions of diquark 
and antidiquark. In the option B in which the a(980) is a KK molecule a Gaussian form for the motion of these 
particles is not appropriate. In this case we take 

i>a(r) ~ - exp (-fir), fi = ^E b m K , (2.25) 
r 

where Eb stands for the "binding energy" of the system. 2 This radial wave function is valid at large distanances 
where the interaction between the two mesons can be neglected. To take it also for small distances is certainly an 
oversimplification. It will nevertheless give us a qualitatively correct picture since at small distances of the center of 
masses of the mesons the constituent quarks are still distributed over the range of w 1 fm. 
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By using Gaussian wave functions to describe the internal structure of the clusters inside the a-meson (i.e. diquarks 
in the genuine tetraquark option and kaons in the molecular option) we find explicitly the form factor g a +^„ 3 (Q ) by 



1 We notice that relati vist ic quark models do not face such a problem. They provide a good description simultaneously of form factors 
and decay constants llSl . 

2 The binding energy of a bound state built up of several constituents is the difference between its mass and the sum of the constituent 
masses. For a bound state in a two-channel problem (e.g. the K°K° and K + K~ channels) one cannot define a binding energy since 
the constituent masses are different in different channels. Nevertheless, we can speak also in this case about "binding energy" through 
the relation with the fall-off of the wave function at large values of r. 



6 



integrating Eq. (| A. 1 2|> 

9a^ v M 2 ) =gexp(-q 2 a 2 D /4), (2.26) 

where g is a ^-independent constant and is the diquark/kaon size parameter. Thus, the ^-dependence of the 
form factor g^q 2 ) is fully determined by a single parameter - the diquark/kaon size. The quantity g — g^ — 0), on 
the other hand, is a function of all the size parameters a a , an, and a rig . 

Strictly speaking, for the decay a(980) — > nr\ a relativistic treatment is necessary because the velocity of the 
outgoing r\ meson is not small. Thus, our nonrelativistic calculation for the a(980) — > 7r?7 decay is not precise, but 
still qualitatively acceptable. 

Numerically, we find for the amplitude of the isovector 1 = 1 a-meson 

A(a+ -► f]n + ) = 3. g sin0 ( m u + m o \ MeVj (2 2?) 



M a - m v 

where is the ?7-meson mixing angle, see Appendix A for details. Respectively, the partial width reads 

r(a° -> ryTr ) - r(a+ -> m +) = 54 I m u + m A ' g 2 MeV _ (2 28) 

\ M a - rrin J 

In the value T(a° — > 7771"°) = 60 ± 13 MeV was obtained making use of the measurements of the full a- width ^| 
and the branching ratio quoted by Particle Data Group [2p| |. 

Our result for the partial width now depends on the values of the constituent quark masses: If we make use of the 
relation (|2.8p . mp + mo = M a — m„, which gives vtiu = to^ = 220 MeV, then g ~ 1 is needed to be compatible 
with ~ 60 MeV for the width. For the values of the constituent quark mass mjj = mo = 330 MeV, sometimes used 
in nonrelativistic quark models, one would need g ~ 0.65. Fig^ exhibits the form factor g for the two scenarios 
as functions of the tetraquark size parameter d - the root mean square distance between the centers of mass of the 
diquarks/if -mesons in the tetraquark. For the Gaussian wave function with the parameter a, one obtains d = -\/3/2 a; 
for the molecular wave function l|2.25|l one finds d = l/(>/2[i). We note that the form factor g does not depend on 
the quark masses. 

In scenario A, the magnitude of the form factor depends strongly on the average separation of the diquarks. Since a 
full overlap of the diquarks would require an unphysically large binding energy we need only to consider the behaviour 
of the amplitude to the right of the maximum. To have g ~ 0.65 4- 1 requires therefore a relatively large distance 
between the diquarks of about ~ 1 1.5 fm. Since the diquarks are extented objects themself this implies that only 
a tetraquark of a large size can explain the width. 

For the KK-mo\ecn\e scenario the dashed curve in Fig. ^ applies. Constituent quarks satisfying l|2.8|) lead to a 
width below 30 MeV for any value of the molecule size. The value for the a — > rjTT decay rate of 60 MeV can only 
be obtained for constituent quark masses around 330 MeV, and requires a molecule with a size of about 0.5 -j- 1 fm. 
This corresponds to an equivalent "binding energy" Eb (d = \j \J2E mK) in the range 40 < E < 150 MeV. The mass 
values M a a = 985.1 ± 1.3 MeV H3, M K+K - = 987.4 MeV, and M ROR o = 995.2 MeV show that an interpretation of 
the a as a K°K° molecule of mixed isospin has the largest but still too small binding energy. In this case, however, 
the r given by (|2.28|l has to be reduced by a factor 2, leading to the values incompatible with the observed width 
of rj 60 MeV. This makes the molecular interpretation of the a° unlikely. A measurement of the decay of the a + is 
needed to shed more light on this question. On the other hand, on the basis of our results, we cannot exclude an 
a(980) structure, in which the K mesons form a molecule at the surface region only, while the interior has a different, 
perhaps two-quark, composition (cf. |2l|). 

2. a=> KK 

For this reaction we consider the decay of the aP particle within the diquark option only. It proceeds via rearrange- 
ments of the constituent quarks. The decay of a -ftT if -molecule into two kaons would require knowledge about the 
formation process of this particle and will not be treated here. 

Within the genuine tetraquark scenario, it turns out to a good accuracy that 

9a^ K {q 2 ) - 0.9g a+ -,, h {q 2 ). (2.29) 

However, different values of q 2 have to be applied in the different reactions. For the a —> KK amplitude, we set the 
momentum transfer equal to zero and find 

K°K+ ) = 3 9 n ( 

_ M a - m K 



A(a+ -> K°K+) = 3.9.9 ( + ™ S ) MeV. (2.30) 
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Fig. 1: The a — > rj-K transition: The form factor g defined by 12.261 for P — rj s vs the size d of the a(980) for the two scenarios: 
genuine tetraquark (solid line) and molecular (dashed line). 



For the ratio of the amplitudes, this gives 



A{a+ V7T+) ~ 0.9 sintf ( 2E±I2£] s 0.72 sintf, (2.31) 



since the ratio of the quark masses (mjj + m£))/{nijj + ms) ~ 1/1.25 is weakly sensitive to their specific values. 

The a-meson is below the K K threshold, and the decay a — * KK proceeds through the finite a-width. Therefore, 
the determination of the branching ratio a — * KK is involved: one must fit the data making use of the coupled-channel 
formula 0. 



3. THE EMISSION OF VECTOR MESONS AND THE VECTOR CURRENT 



A. The vector current for constituent quarks 

A proper vector current can interpolate vector meson fields. In the example of an isovector current one has for 
the p meson field operator <I>^ = — uj^d. As in the case of the pseudoscalar current, the current for constituent 
quarks should no more contain the p-meson pole which occurs in the current formed by the current quarks of QCD. 
In the q 2 region around and below the p meson resonance one can write 

9 

Tfl — 

u^d=—^U^D. (3.1) 
□ + trip 

The normalization of the right hand side is simpler than in the pseudoscalar meson case since it must be fixed to be 
one at momentum transfer q 2 = 0. Thus, the interpolating p-mcson field becomes 

TY) 1 — 

= T^pr— - U^D. (3.2) 
p fv U + mf, 

The amplitude for a decay proceeding by the emission of a p~ meson can thus be obtained by sandwiching the 
constituent quark operator 



T=-^e;(q)UrD (3.3) 
between hadrons formed by constituent quarks. Considering now the vector analogue of l|2.5|l 



U^D = iU T ~f°j 2 ^D, (3.4) 
UD = iU T f^ 2 D, 
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one gets for conventional transitions in nonrelativistic approximation 

jo = u*D, (3.5) 

Jv = ^^—{U^acJkdkD - dkU^kSD). 
zniQi 

For the fall-apart operators one finds on the other hand 

J v = ~7^— d k {U T a 2 cT k D), (3.6) 
Jv = -i(U T a 2 aD). 

It is seen, that the current operators for fall-apart transitions are particularly simple. The corresponding matrix 
elements can be expressed by overlap integrals without derivatives. The divergence of J M has to vanish and gives 
immediately the constraint for the effective quark masses 

Mi - M f ~2m Q . (3.7) 



B. Fall-apart amplitude for spin-1 polyquark mesons 

We consider the fall-apart process X(l + ) — * V(l~) with the emission of a vector meson. The case in point here are 
the decays A"°(3872) — > J/ip -kit and X°(3872) — > J/ip iririr, mediated by the p° and ui meson, respectively. Clearly, 
the isovector component of X a contributes to the first reaction, while the isoscalar component contributes to the 
second one. 

Let us briefly outline the procedure for the isovector X° transition. We start with the meson transition amplitude 
induced by the conserved vector current j v — -^(uj^u — dj^d), and write its decomposition as follows 

(V(p')\jV\X(p)) = e^hiq 2 ) + (p„ ■ q 2 - % ■ qp)e pp , ££ , f 2 (q 2 ) + (ep')W e> hrf) + (e'p) e wpl£ fi{q 2 ). (3.8) 

The main contribution to the decay rate of the reaction X°(3872) — * J/ipinr comes from the region where the 
intermediate /o°-meson is nearly on-shell. In the X-rest frame, the intermediate p° is almost at rest. Therefore, we 
can neglect the form factors f 2 , /3, and f 4 in Eq. (|3.8|l and keep only the form factor f\. The form factor fx contains 
the p°-pole so we may write 

o 

711: 

h(q 2 )= ' ^ (q 2 ). (3.9) 

— q A + m z p 

The amplitude of the X° — » J/ip p° transition then takes the form 

A(X a -> J/rP p°) = ^e;(q)e:(p')e,(p) fo ^Fi(m»). (3.10) 
jv 

Making use of the relation (|3.1|l . the form factor F\ may be obtained from the amplitude of the constituent-quark 
vector current 

(V(p')\^(U^U-D^D)\X(p))=e^e* u (p')e x (p) q aFl {q 2 ) + ■ ■ ■ , (3.11) 

where • • • denote small terms containing higher powers of the small momentum q. The z-component of this equation 
is sufficient for calculating Fi(q 2 ): 

F x = ±(Vtf = -q^^iU^U - D^D^Xip^ 0,±)). (3.12) 

<r V2 

The ± signs in the state vectors refer to the particle polarisations. In the X-rest frame, the J ftp is moving slow, and 
a nonrelativistic approach may be used for the calculation of the form factor. Further details of the calculation are 
given in Appendix A. 
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C. Fall-apart decay of the axial- vector polyquark X(3872) 

The recently observed charmonium-like V(3872) particle is likely a J PC = 1 ++ state. Its mass Mx = 3871.3 ± 
0.7 ± 0.4 7], its small width T(X) < 2.3 MeV and its decay properties make it a good candidate for a polyquark 
hadron. Like in the case of the a particle we consider two options for this state: 

A: X is a confined tetraquark consisting of two color-triplet diquarks in a relative S'-state, one diquark (anti-diquark) 
with spin and the anti-diquark (diquark) with spin 1. 

X q = (Q T io 2 C) (Q T io 2 a C) + (Q T io 2 o C) (Q T ia 2 C). (3.13) 

B: X is a four-quark molecular state: a weakly bound 5-state of a pseudoscalar D meson and a vector D* meson 

X q = (C T ig 2 Q) (C T ia 2 a Q) + {C T ia 2 a Q) (C T ia 2 Q). (3.14) 

For the transition A°(3872) — > J/ip p° mainly the isovector component X I=1 — -^(X u — Xd) contributes, whereas 

for the tranition V°(3872) — > Jji\> oj it is the isoscalar component X I=Q = ^(X u + Xd). The physical X will in 

general be a combination of X I=1 and X I=0 . The Fock state representations of the X and of the J/ip are given in 
Appendix A. For the wave function of the confined diquark- antidiquark system we take a Gaussian form 

$ x (r)~exp(-r 2 /2a|), (3.15) 

with r the distance between the center of masses of the diquarks. 

In case B, we have a bound state of two colorless objects. Since the properties of a weakly-bound state are largely 
determined by its binding energy Ef,, we take for the relative motion of the two constituents the wave function 

$x(r) - - exp (-nr), fi = y/2E b M D M D ./(M D + M D .). (3.16) 
r 

The X meson mass is close to the D°D*° threshold at 3871.6 MeV, and around 7 MeV below the D+D*~ threshold 
at 3879.4 MeV). Thus, the X binding energy is restricted to the range E b < 7.5 MeV. 

It is convenient to express the decay rates via the transition form factor (using nonrelativistic normalization of the 

states) 9x^j H} = Fx ™*- j ^ x . 



9x^{f) - (J/V =1 ' J ' =1 (-Q)\^(0 T ^2<7 3 U) - (D T ia 2 a 3 D))\X J = l ' J ' =1 (0)). (3.17) 

It turns out that the result for gx->j/ip expressed in terms of the radial wave functions of the composites, for the 
diquark-antidiquark option is \/3 times larger than the analogue result for the molecular option. In other words, if 
the spatial distributions would be the same in the two options, the widths would differ by the colour factor 3. 

For numerical estimates we use the following inputs: The U C scalar and vector diquarks, as well as the D and 
D* mesons are described by Gaussian wave functions with the same size parameter ajj = 0.6 fm. For J/tp we take 
a J/ip = 0.5 fm and set Mx = 3872 MeV and Mjm = 3097 MeV. The values of other relevant parameters, as well as 
the equations related to the finite widths of the p and u>, are given in Appendix A. 

With Gaussian wave functions for the structure of the J7C-diquarks/-D-mesons forming the polyquark X, the 
transition form factor has the form 

9x^j/^{q 2 ) = gexp(-g 2 a|,/4), (3.18) 
where ao is the UC diquark/D-meson size parameter. Then the rates obtained are 

r(V /=1 -> J/^tt+tt-) = 5.2 (^) 2 MeV, 

T(X I=0 -> J/^jn+TT-TT ) = 1.4 (t^) 2 MeV. (3.19) 

Taking into account that the branching ratios of the two decay modes seem to be close to each other 

Br(X -> J/i/;7r + 7r-7r ) 



Br(X -> J/^tt+tt-) 



1.0 ±0.4 ±0.3, (3.20) 
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Fig. 2: The X — > J/ip transition: The form factor g defined in 13. 1811 vs the X size parameter d for the genuine tetraquark 
(solid line) and the molecule (dashed line). 



we conclude that X(3872) should be dominantly an isosinglet particle 

X = cose x X I=0 + sin9xX I=1 , sin 9 X ^ 0.46 ± 0.3. (3.21) 
For the central value of the mixing angle 9x we find 

T(X -► J/^m+TT-) ~ T(X -> J/V-Tr+TT^Tr ) = 1.1 (^) 2 MeV. (3.22) 

The dependence of the coupling g on the V-size parameter d - the root mean square distance between the center of 
masses of the two clusters inside X - is plotted in Fig. [21 We only need to consider the regions for d to the right 
of the maximum. Lower values are very unlikely: they would correspond to a strong overlapping of the diquarks/Z? 
mesons and consequently to an unobserved large binding energy. Thus, only a large separation can reasonably well 
explain the small width, which must be below the total width of ~ 2.3 MeV. In the diquark-antidiquark scenario the 
necessary small width can be obtained for an averadge distance of the two clusters equal or larger than d > 2.5 fm. 

A pure molecular picture would require an even much larger size corresponding to an extremely small effective 
binding energy Ef, < 0.2 MeV (related to d according to d = \/\f2Mj}Eb). Such a state is not likely, but still 
conceivable in principle, since the sum of the masses of D° and D*° coincides with the mass of the X within error 
limits. We can exclude X mesons composed purely of charged D mesons (D + D*~ + D~D* + ). The corresponding 
binding energy obtained from the mass values is ~ 5 MeV which would lead to a much too large width T ~ 10 MeV. 

We therefore conclude that the diquark picture is preferred. In any case, a polyquark hadron must possess an 
unusually large spatial extention in order to have a small decay width. 



4. FALL-APART DECAYS OF POLYQUARK BARYONS 



The observation of a pentaquark at 1530 MeV could not be confirmed. However, exotic baryons of higher mass and 
pentaquarks containing heavy quarks could still exist. For a decay into a conventional baryon under the emission of 
a pseudoscalar or a vector meson a treatment analogous to the one given in Sections 2 and 3 can be performed. The 
calculations for the pentaquark 8(1530) given in Refs. 0,H3,|3| are n °t repeated here. It was found, that a small 
width of the order of 1 MeV requires a large spatial extention (molecule like) of the pentaquark, similar to what we 
found here for polyquark mesons. We collect in appendix B formulas for pentaquark decay amplitudes which are more 
detailed than the ones contained in the mentioned references. These equations may become applicable in the future. 
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5. CONCLUSION 

We studied fall-apart decays of polyquark (tetra, penta) hadrons within the constituent quark picture. By making 
use of chiral symmetry the axial current for constituent quarks is shown to contain a local part Qj^j^Q and a nonlocal 
contribution proportional to the pseudoscalar density Q75Q. Using also the close connection between vector currents 
and vector mesons it turns out that fall-apart processes with the emission of pseudoscalar mesons and vector mesons 
can be calculated from simple overlap matrix elements. The transition amplitudes depend then little on the relative 
velocities of the constituent quarks, but rather decisively on the size of their spatial distribution. These facts can help 
to study the states suspected to be polyquarks and the problem of the clusters they are made of. 

We applied the developed formalism to the analysis of the decays a° — ► i]tt, KK and X — > J/ipn + ir~ , J / 'ipir + tt~ ir° , 
assuming that a and X have a polyquark structure. We tested two extreme scenarios, namely: (A) a and X arc 
confined diquark-antidiquark states (genuine tetraquarks) and (B) these states are bound states of two i-T-mesons and 
two D-mesons, respectively, i.e. molecule like particles. We calculated the decay rate for both scenarios as a function 
of the averadge distance between the center of masses of the two clusters. 

I. For the a (980) we found: 

• Within the scenario (A) one can reproduce the transition rates correctly, if the averadge distance between the 
two diquarks is taken to be relatively large, around 1 -e- 1.5 fm. Since both diquarks are extented objects themself, 
the a will be a relatively large-size tetraquark object. 

• The situation for the a , if described as a KK molecule, is more complex: If the two if-mesons essentially keep 
their identity in the bound state, the aP would have to be of the isospin mixed form K°K°. The calculated width 
for the process a — > ijtt then disagrees with the observed decay rate making this picture for the a unlikely. On 
the other hand, if the KK structure refers only to the outer part of the particle, and the molecule description 
is not valid for the inner region, the experimental transition rate can be accomodated for the average distance 
between the clusters around 0.5 -j- 1 fm. The behavior of the wave function in the outer region would then 
correspond to a (fictitious) "binding energy" larger than 40 MeV. 

II. For the X°(3872) meson we found: 

• Within the scenario (A) (X consisting of two charmed diquarks) we managed to obtain the desired small decay 
widths which are compatible with the experimental limits. Necessarily, the X(3872) must be mainly an isoscalar 
1 = particle: 

X = cos6 x X I=a + sin6 x X I= \ sin 6 X ^ 0.46 ± 0.3. (5.1) 

This conclusion is based on our estimate for the X 1=1 — ► J/tp tttt and X I=0 — ► J/tp tttttt decay rates and the 
measured ratio of the T(X — > J/tp tttt)/T(X — ► J/tp tttttt), which has so far a rather large error. To obtain 
sufficiently small widths the X° must be a large size particle. The distance parameter d has to be around or 
larger than 2.5 fm. 

• For the DD*-molecule scenario our finding is again somewhat involved. In case the D mesons keep essentially 
their identity in the bound state, this X meson would be a particle of mixed isospin (i.e. made of zero charged D 
mesons only). Then, the binding energy would have to be very small, smaller than about 0.2 MeV corresponding 
to a huge radius. This makes this case unlikely, but still conceivable in principle, since the sum of the masses of 
D° and D*° coincides with the mass of the X within error limits. 

Evidently, if polyquark hadrons exist, these particles should be rather large-size objects. A similar result was obtained 
in the earlier discussion of the exotic baryon where we also found that a small width for fall-apart processes is correlated 
with a large particle size. 

Acknowledgments. We thank V. V. Anisovich, P. F. Ermolov, D. Gromes, W. Lucha, O. Nachtmann, and S. Simula 
for useful discussions on this subject. D.M. is grateful for hospitality and financial support from the Institute of 
Theoretical Physics of the Heidelberg University during his stay in Heidelberg. D.M. was supported by the Austrian 
Science Fund (FWF) under project P17692. 
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APPENDIX A: NONRELATIVISTIC FOCK STATES AND TRANSITION AMPLITUDES 

Note: In this appendix A all quark fields refer to constituent quarks. For a convenient and lucid presentation we 
will use here small letters to denote them (in contrast to our notation in the main part of the paper where we had to 
distinguish between current and constituent quarks). 

1. Decay of scalar tetraquarks to two pseudoscalars 

a. The light mesons 

The rj and 77' mesons are mixtures of strange and nonstrange components |2fil | 

|t?) = cos%„) - sin6>|?7 s }, 

\rf) =sin%„) +cos% s ), (A.l) 

where r\ n = (uu + dd)/y/2, rj s = ss and sin# = —0.65. For the amplitude (ri\d T ia 2 u\a + ) only the strange component 
of the ?y-meson contributes. Thus one has (r]\d T ia 2 u\a + ) = — sin# (il s \(d T i<7 2 u)\a + ). 
The r) s component has the structure (summation over colour is implied) 

, ... 1 f ( Tti\T\ + m 2 r 2 \ _ . , . , T . 

(r] s (p)\ = —= I dridr 2 exp \-ip ■ ) <P Vs (n - r 2 ) \(s 1 (n) ia 2 s(r 2 ) . (A.2) 

V6 J \ mi+m 2 J 

The radial wave function $ rh (r) is normalized according to 

cF|$ r)s (r)| 2 = 1. (A.3) 

We shall use the Gaussian parametrization $(r) ~ exp f — go 2- )' ^ or ^ ne ^ s s * a t e one nas to take mi = m2 = m s . 

The wave function of the if -meson has the same form. Here one has to put mi = m s and m 2 — and to replace 
a r)s by ax- For the pion one should set mi = m 2 = m>d and a Va — > a n . 

b. The scalar tetraquark a + 

We consider at first the scalar tetraquark meson a + as consisting of a spin-zero diquark us and an anti-diquark sd 
in a relative L = angular momentum state. This tetraquark Fock state has then the following form 

1 + / -»\\ /5 / j- j- j- {.-.msn +mr 2 + m s r 3 + mr 4 \ _> 

a W) = v3 / dridr 2 dr 3 dr±exp ip — ■ W Q + (ri, r 2 r 3 , r 4 ) 

J \ 2(m s + mj J 

xD^in^DU^nm. (A.4) 

.D a denotes the bilocal diquark annihilation operator 

D a (n,f 2 ) = 4|e aaiQ2 (s aiT (ri) «t 2 q a H?2)) , (A.5) 



and I? a the corresponding anti-diquark operator. The diquark picture requires the coordinate wave function of the 
tetraquark to have the factorized form 

V a +(ri,r 2 \r 3 ,r 4 ) = ^D(ri 2 )^ D (f 3 i)^a + (p), (A.6) 

where 

m s ri + m d f 2 



r\ 2 = f\ - f 2 , R\ 2 = 
^34 = f 3 - r 4 , i? 34 = 



m 3 + m d 
m s r 3 + mdVi 



m s + m d 

p = R 12 - i? 34 , (A.7) 
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and 



J |$ D (r)| 2 dr = 1, J |$ a+ (p)| 2 dp= 1. (A. 



Again a Gaussian parameterization for the wave functions is used: <j>.o(r) ~ exp ^ — 2a 2 ") anc ^ 0^) ~ ex P ( — 2a 3- ) ■ 
As an alternative, we take a + to be a KK molecule. In this case, (|A.4J| has to be replaced by 

i + / -\\ /"j-j-j-j- (■ - +mf 2 +m s r 3 + mr 4 \ 

a W) = / dridr 2 dr 3 dr 4: exp ip — ■ W a + (r 1( r 2 r 3 , r 4 J 

xif+t(f=i,r' 2 )^ t(r'3 ) f 4 )|0) ) (A.9) 

with 

^ + (ri,r 2 ) = i (s T (ri) «CT2 u(rb)) , 

K°(r 3 ,n) = ±= (dT(r 3 ) ia 2 s(f 2 )) . (A.10) 



c. The a + — > 7771- + and a + — > K + K decays 

We introduce the dimcnsionlcss form factor g a +^p (in nonrclativistic normalization for hadron states) where P 
stands for rj s or if : 

Sa+^p(<T 2 ) = g Q A (P(~q)\(f i°2 u)\a + (p= 0)). (A.ll) 
Simple algebra leads to the relation 

' TOir'i — TO 2 r' 2 



g a +^p(q 2 ) = -9a k / dridr 2 exp (zg miri +m2?2 j $ p (^ _ ^ q> D (f 1 ) $ D (r 2 ) $ Q 
7 V m 1 +m 2 J 



mi + m 2 



(A.12) 



The indices 1 and 2 correspond to the constituents of the final meson. For r\ in the final state these are s and s, for 
the K° meson s and d. The value of n is in case the a + is a diquark composite while one has k = -4g if it is a 
KK molecule. 

The transition amlitude for the process a + — ► ?y7r + is then given by 

, 2 ,m v + m D 



A(a + — > mr + ) = sinfl 5 a +_ t7J<> (g ) u ; " vMM^. (A. 1 3) 

For the a —> r) transition, the equal quark masses drop out from the form factor g Eq. I|A.12|) . So the only place 
where the quark masses come into the game is the factor mjj + mr> in the amplitude. 
Neglecting the a + width, the decay rate a + — > 7]7r + takes the standard form 

r(a+ -> T77T+) = — I J ^ 1 smfl g a +^r ls (Q )| • (A.14) 

In the above formulas -E 1 ^ and — <f are the energy and the spatial momentum of r\ in the o + rest frame; and m,r 
denote the masses of rj and tt, respectively. In this process the ir field is used as an interpolating field. 
Similarly, the amplitude of the a — > KK decay reads 

A(a+ K + K a ) = g a ^ K (<f 2 ) mU + mS y/4M^. (A.15) 

JK 

According to our findings, to a good accuracy g a +—,K(q 2 ) — 9a+—>r)s(<l 2 )- 



14 



2. Decay of the axial- vector tetraquark to two vector mesons 

a. The "genuine" tetraquark X q 
For the tetraquark state, we take the following nonrelativistic representation: 
-;=i ; =i/.^\ v 7 ^ f ,_, m c ri + mr 2 + m c r 3 + mr 4 



^ = ^ dridr 2 dr 3 dr 4 exp[ip^^ L . ] $x(Ri3 - R&) 

V 2 J \ 2(m c + mj y 

x (D?^ (n, r 2 )D$(r 3 , r 4 )$ Dl (^2)^(^4) + ^(ri, rsjsftfo, u)$ Do (r 12 )<S> Dl (f 34 )|0) 

(A.16) 

where the bilocal diquark annihilation operators are defined as follows: 

Dg(n,r 2 ) = 4^ e<laia2 (« aiT (^i) ^2 c Q2 (r 2 )), 

v 12 



Oi,i(r-3,r 4 ) = ^ eQa3a4 (<? Q3T (^ 3 ) ^25 fa - ia 2 ) c a *(f 4 )). (A.17) 

I?0j -Di denote the corresponding anti-diquark annihilation operators. According to the diquark picture of the 
tetraquark given in l|A.16|l . the factorized tetraquark wave function contains $d and $_Di, the radial wave func- 
tions of the spin-0 and spin-1 diquarks. $x stands for the wave function of the confined bound state composed of the 
two diquarks. 

b. The DD* molecular state X q 
For the molecular state of the X q the following nonrelativistic representation is taken 

\X J q - 1 '^-\p)) = — / dndr 2 drzdr A exp Up \, \ \ \§x (R12 - #34 

y V2 7 V 2(m c + m) J 

x (Ct!i(^i^2)^ t (r 3 ,r 4 )$i )l (ri 2 )$i )0 (r 3 4) + ^(r!, r 2 )^a(^3, r 4 )$D (^12)$^ (f 34 )) |0) 

(A.18) 

where 

D(n, r 2 ) = -±=<5 Q1Q2 (c° ir (ri) l( r 2 r (r 2 )) , 

£> 1 * 1 (r 3 ,r 4 ) = -L,5 Q3 ^(c a3T (r 3 ) ia^fa -ia 2 ) <T(r 4 )), (A.19) 

$£> and are the radial wave functions of the spin-0 and spin-1 mesons and $x is the molecular wave function 
of the bound state composed of the two mesons. 

c. The J/tp state 



(J/iP J=1 ' J > =1 m = % / dndf 2 exp Lip^±^j ^(fi - r 2 ) (0|(c a < T (ri) ^fa - ia 2 )c a '(r 2 )). (A.20) 

d. The X — » J/tp transition amplitude and the X — > J/ip-Kix and X — > J/iPtttyty decay rates 

Let us consider the fall-apart process Af(l + ) — > Af(l~) with the emission of a vector meson. A case in point is the 
decay A°(3872) J/ijj p° -> J/^> tt+tt" and AT°(3872) -> J/V> w -> J/^> tt+tt'tt . 
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To obtain the transition amplitude of the isovector component of A (3872), which we denote A /_1 , we start with 
X 1=1 — > J/ip transition induced by the conserved isovector vector current 

j v = ^=(u^u-d^d). (A.21) 

The form factor decomposition reads 

(V(p')\j%\X(p)) = e M ee'h(q 2 ) + iPn ■ q 2 - ■ qp)e pp , ee ,f 2 {q 2 ) + {ep')t PPP >s> h{q 2 ) + {e'p)e WP 'eh{q 2 )- 

(A.22) 

The main contribution to the decay rate of the reaction A (3872) — > J/ip tt + tt~ comes from the region where the 
intermediate p°-meson is nearly on-shell. In the A-rest frame, the on-shell p° meson is produced almost at rest. 
Therefore, we can neglect the form factors f 2 , /3, and f± in Eq. (|A.22|I and keep only the form factor fx. The form 
factor fi contains the p°-pole so we may write 

2 

171 

hiq 2 )-^f^F,{q 2 ). (A.23) 



The amplitude of the A 7-1 — > J/ip p° transition then takes the form 

mp 
fv 



A(X I=1 - J/4, p°) = ^e;(q)e:(p')e x (p) q^^F^m*^ (A.24) 



with fy defined by 



(0|j>°) = e p f v m p , f v = 216 MeV. (A.25) 
Treating the p° as a stable particle, one finds for the A — » J/ip p° decay rate 

T(X'=^ P °JM = - ^J xf f (A.26) 

where q is the momentum of the p meson in the A-rest frame, q 2 <C m 2 . However, the p° is unstable leading to 
the final state ir + ir~. Taking into account the finite width of the p°-meson, we obtain for the A /=1 — * J/tp tt + tt~ 
amplitude 

A(X I=1 -» J/V^+tt-) = A(X I=1 -> J/^P°)— 9 l ———A{p -» tttt), (A.27) 

m 2 p - s- B pp 

where B pp {s) is the p-meson self-energy function, explicit expression for which is given in |25|. 
The p — > tt + tt~ amplitude may be parametrized as 

A(p -> tttt) = -g p ^{k - fc%e", (A.28) 

e M is the p- meson polarization vector. The corresponding decay rate for the virtual p-meson with the mass i/s reads 

q 2 ( 4m 2 ^ 3/2 



and 

lmB pp (s) = ^T p (s). (A.30) 
The decay rate of the reaction A — > J/i/jtt + tt~ takes the form 

dT(X I=1 -» J/^tt+tt-) 1 fl 2|jr 1 ( fl )|2A 1 /2(M=,M2, a ) ^r p (s) 



(A.31) 
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where ReB pp (s) can be found in [23. Let us notice that if we take the limit T p (s) — > (i.e. g p7r7r — > 0), the decay 
rates satisfy the simple relation T(X I=1 — > J/-07r + 7r~) = r(A /=1 — * J/ipp°). For numerical estimates we use the 
values m p o = 773.8 MeV and g p7r7r = 11.4 from the recent analysis j2^|. 

We now calculate the form factor F\ within the constituent quark picture. Making use of the relation (|3.1|) . the 
form factor F\ can be obtained from the constituent-quark vector current 

(V(p')\-^=(tJ^U - Lhr M D)\X(p)) = e^et(p')e x (p) q^tf) + ■ ■■ , (A.32) 

where • • • denote small terms containing higher powers of the small momentum q. The z-component of this equation 
is sufficient for calculating F\(q 2 ): 

Fi = ^(Vtf = -q^^iU^U - D 1 3 D)\X(p = 0,±)). (A.33) 



The ± signs in the state vectors refer to the particle polarisations. In the X-rest frame, the J/ifj is moving slow, and 
a nonrelativistic approach may be used for the calculation of the form factor. 



Isolating the kinematical factor related to the normalization of hadron states, we can express F\ by gx^j lip{q 2 ) 



W) = 9x^j/4,{q) I f „_;,, , (A.34) 

with 



M x - E JM 



9x^j,4f) = (JH^ 1M=1 (~q)\-^{{U T 1*2*3 U) ~ (D T ta 2 a 3 D))\X J=l ' J ' =1 (0)) • (A.35) 
Here the standard nonrelativistic normalization of states is used. Explicit calculations lead to the expression 

gx^.;Mq 2 ) = (J/V J=1,J " =1 (-«)I(Q ' T (*'= 0)0-1 Q(r = 0))|X J = 1 ^= 1 (0)) 

= -K J dr^exp fiq rx + r2 ^ ^(rf )$ Do (ff )$ x ( — (^1 - r 2 )\ Qj/^n - r 2 ), 

(A.36) 

where k = 1 if X is a diquark composite, and n = -^j if X is a molecule formed by D and D* mesons. Here q is the 
momentum of the outgoing J/tf) in the X-iest frame. The form factor gx^j/ip 

(q 2 ) determines r(A /=1 -» J/^+O 

decay. 

A similar treatment is applied to calculate the three-pion decay X°(3872) — * J/i/j tttttt via the w meson. In this case 
the isoscalar component X I=0 determines the amplitude. The corresponding width T(X I=0 ~ + J/ijj tttttt) is obtained 
using the same form factor gx~>j/^(q 2 ) by a formula similar to (|A.31|> with r p (s) — * T w (s) and m p — * m u , and 
multiplying by the Br(cu — > 3tt) = 0.89. Because of the small width of the w-meson, the s-dependence of r w (s) makes 
little difference, mainly the value r w (m^) = 8.5 MeV is essential. 



APPENDIX B: FALL- APART AMPLITUDES FOR SPIN-1/2 POLYQUARK BARYONS 

In this section we discuss the baryon to baryon transition matrix elements induced by the axial-vector current and 
the corresponding decay amplitudes for the emission of a light pseudoscalar meson. 

As an example, we consider the — > N transition amplitude induced by the strangeness-changing axial current 
S7/i75U, where is an exotic polyquark hadron and N denotes a conventional spin 1/2 baryon. The corresponding 
hadronic decay is — ► NK, The amplitude of interest has the following general decomposition in terms of invariant 
form factors 

(N(p')\s'y l SY5u\Q(p)) = g A {q 2 )u N (jp')-i^ b u e (jp) + gp{q 2 )q ll u N {p')^u e {p) + g T (q 2 )u N (p')(J^q u 'l5Ue(p) (B.l) 

with q = p — p'. Since the K pole of interest occurs in the form factor gp{q 2 ) we can define a residuum function r(q 2 ) 
by setting 
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Taking the divergence of the axial vector current as an interpolating field for the K-meson, the decay amlitude is then 
given by 

T(0 -> NK) = geNK • u N {p')ilbUe(p), (B.3) 

with 

geNK = r -^M. (B.4) 

IK 

In the chiral limit with ttik — > the axial vector current is conserved leading to a relation between g A (q 2 ) and r(q 2 ) 

(Me+Mjv) g A (q 2 )=r{q 2 ). (B.5) 

One expects that gA and r do not change significantly by going to the chiral symmetry limit. Therefore, we have two 
possibilities to calculate gsNK in a model for hadrons, namely from r(ro^) and from <ju(m^-) |10| . 

By expressing the axial current by virtue of eq. 1)2. 1[) in terms of the constituent quark field operators one gets 

(N(p')\sr l5 u\e( P )) = g^(N(p')\S^ l5 U\e(p)) - (ms + ma) _ X 2 g2(N(p')\S 75 U\Q(p)). (B.6) 

The matrix elements on the right-hand side can be expressed in terms of invariants in the same way as done above. 

(N(p')\Sn 5 U\e(p)) = G A {q 2 )u N {p')Y l l 5 u B (p) +G P (q 2 ) q »u N (pi) l5 u e (p) + G T {q 2 )u N {p')o^q vlb u e (p), 
(N(p')\S l5 U\Q{p)) = G 5 ( q 2 )u N (p') l5 u e (p). (B.7) 

Now, however, all form factors are regular functions and have no poles in q 2 in the g 2 -region of interest. 
By comparing l|B.7(l with (|B.1|) we find that the form factors are related to each other as follows 

9A(q 2 ) = 9 Q A G A {q 2 ), 
g T (q 2 ) = g%G T (q 2 ), 
r(q 2 ) 



- q * + m% =«^ )_ =?+^ W) J- (R8) 
At the pole one has 

r{m 2 K ) - g Q A {m s + m v ) G 5 (m 2 K ). (B.9) 

The divergence equation for the first part of the axial-vector current of constituent quarks l|2.2|l leads to the following 
relation between the form factors 

(Me + M N ) G A (q 2 ) - q 2 G P {q 2 ) = (m s + m C /)G 5 (g 2 ). (B.10) 

This relation is automatically satisfied in the relativistic dispersion approach of Ref. |27|. In general however, 
calculated with trial wave functions for initial and final hadrons, (|B.10|) is not automatically satisfied. 
At q 2 = m 2 K we have 

r(m 2 K ) = g Q A ((M e + M N ) G A (m 2 K ) - m 2 K G P {m 2 K )) . (B.ll) 



The form factors Gi can be calculated from different components of the l.h.s. of 1B.7|I for different polarizations of 
the initial 0. We work in the rest frame p = (Me, 0) and choose q = (q^, 0, 0, |q|). It is convenient to use now the 
nonrelativistic normalization of the state vectors. Then the form factors are given by the equations 

(M N - M&)(E N + M N ) Mjv + M© M| + M% — 2M®E N 

A — A H At. H ; : [Ar — At ), 

A 2M e \q\ 2M@ L (E N -M N )2M @ v h 

\q\ E N - M N 

Gs _ (Ej^MA^ (B12) 
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with 



a° = (i\rV)|S7°75tf|e T >, 

A T = i(iV T (p')|57 2 75{/|e i ), 

A 5 = (N^jf)\S l5 U\^), (B.13) 



p 1 = —qlies in the negative z-direction, \q\ = ^/-E^ — M%, E N = jM^i-^e + ~ <7 2 ) w i tn <? 2 = m if f° r tne decay 
process. The relation Al = At for q = guarantees that the form factors Gi are finite at q = 0. 

For a transition in which the final baryon moves nonrelativistically we start by writing the form factor decomposition 
appropriate for nonrelativistic motion in the O rest frame: 

(N(-q)\J°\e) = F (&ff?fe), (B.14) 
(N(-q)\r A \Q) = Fii&Jfa + MtWZfte), 



where £jv,e are two-component nonrelativistic baryon spinors. This paramterization gives for the amplitudes in ljB.131) 

A" = \q]F (q 2 ), (B.15) 
A L = F 1 (q 2 ) + \q\ 2 F 2 (q 2 ), 
A T = Ftiq 2 ). 

Taking into account the structure of the quark currents for a fall- apart process (|2.7[) . the amplitudes A and A 5 are 
related to 1/toq accuracy 

A = -A 5 . (B.16) 

The amplitudes Al and At involve derivatives of the wave functions and may thus be sensitive to subtle details of 
these wave functions. The amplitude A — — A5, on the other hand, is a simple overlap matrix element. 
To the accuracy of our nonrelativistic approximation, the solution (jB . 1 2|> takes the form 

Me-Mjy M e + M N {M e -M N )\ 

Ga= 2M^ 2MnFo+ 2Mb Fl + 2Me 

2M N 1 M e + M N 

Gt = ■ — Fc\ — Fi H — — 2 A I \r Fo . 

2M e 2M e 2M @ 

G 5 = -2M N F . (B.17) 

Now we can apply the divergence equation (|B . 10|> . Setting q 2 = (Ms — Mn) 2 — j^q 2 , and neglecting again terms of 

order q 2 /M 2 , we see that the terms proportional to F2 drop out from this equation. As in the meson case considered 
above, this equation reduces to a constraint for the constituent quark masses, namely to 

F1 

M e - M N - — =m s + m u . (B.18) 
^0 
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